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We study the effect of quantum and thermal fluctuations as well as 
the mass dependence of the phase structure of QCD at finite tempera- 
ture and density within a dynamical Polyakov-loop-extended quark-meson 
model. The glue dynamics is simulated by the Polyakov-loop potential, 
also including the back-coupling of the matter sector to the glue dynamics. 
In the chiral limit, the chiral phase transition at large chemical potential 
and low temperature splits into two transition branches. For non-vanishing 
pion masses the chiral transition at small chemical potential changes from a 
phase transition to a crossover. We close with a discussion of a systematical 
improvement of the current model towards full QCD. 

PACS numbers: ll.30.Rd, 64.60.ae, 12.38.Aw, ll.lO.Wx 

1. Introduction 

Over the past two decades, much progress has been made experimentally 
and theoretically in our understanding of the phase structure of strongly- 
interacting matter [1]. This has been achieved in a combination of first 
principle continuum computations in QCD, see e.g. [2], lattice simulations, 
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see e.g. O IH [5], as well as computations in low energy effective models 
for QCD, see e.g. [H O El O [JOl HI]. At vanishing chemical potential it 
is established by now that there are two crossover transitions at similar 
temperatures, see e.g. [HOES]. Presently, it is still much under debate 
whether the two transitions coincide over the whole phase diagram or de- 
viate at some chemical potential, leave aside the phase structure at larger 
chemical potential. In lattice simulations progress is hampered by the sign 
problem at finite chemical potential. The sign problem does not affect the 
continuum approaches which, however, have to deal with the increase of rel- 
evant degrees of freedom and fluctuations as well as the strong correlations 
at increasing density. 

In this context it has been argued that the QCD effective models can 
be understood as well-controlled approximations of full first principle func- 
tional approaches to QCD, see [21[T3]. This connection allows to successively 
constrain and finally determine the sensitive model parameters, and hence 
to gain qualitatively on predictive power within these models. In turn, the 
models provide a more direct access and understanding of the mechanisms 
at work for e.g. the chiral and confinement /deconfinement dynamics, see e.g. 
[131 IS]- The Polyakov-extended quark-meson (PQM) model |7j is specifi- 
cally well-adapted to the embedding into first principle approaches to QCD. 
By now it is well-studied also beyond the mean-field level at vanishing |16] 
and finite density \13\ \T7\ \TE[ [T9] . In this model the matter fluctuations are 
included directly within a fully non-perturbative functional renormalization 
group (FRG) approach, for reviews on QCD and QCD-effective models see 
[201 mi [221 [Ml [Ml [25] • 

In Polyakov-extended models (part of) the glue dynamics is included 
with a Polyakov-loop potential in the free energy. This potential is usu- 
ally constrained with the expectation value and thermodynamics of lattice 
Yang-Mills (YM) theory [9l[T0]. Clearly, this does not completely fix the 
potential and the phase structure is sensitive to the chosen potential, see e.g. 
|26| . This necessitates a better determination of the potential. Corrections 
due to matter fluctuations have been computed in [7] , reviewed in [27j , and 
have been confirmed at vanishing temperature with a direct comparison of 
the Polyakov-loop potentials in Yang-Mills theory [M] and QCD [12]. Fur- 
ther work aiming at the determination also at finite chemical potential is in 
progress. In summary, the PQM model provides us with a good approxima- 
tion of low energy QCD at not too large chemical potential, and its direct 
relation to continuum QCD approaches allows us to strengthen this relation 
successively. 

In the present contribution we study the interrelation between chiral 
symmetry breaking, its restoration at large temperatures and/or chemical 
potential and the confinement/deconfinement transition within the dynam- 
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Fig. 1. RG flow for QCD, including gluon, ghost, meson and quark contributions, 
respectively. 

ical two- flavour PQM model. We map out the full phase diagram and 
also provide a conservative estimate for the possible emergence of a criti- 
cal endpoint. Further insights into the mechanisms at work are extracted 
by varying the explicit chiral symmetry breaking parameter, and hence the 
Goldstone-boson (pion) mass. 

2. QCD flows and the matter back-coupling 

The dynamical Polyakov-extended quark-meson model used in the present 
contribution consists of a Yukawa-type model of propagating quarks and 
mesons (pions vr and sigma a). As a low-energy effective model of QCD it is 
formed dynamically when fully taking into account the four-fermi interac- 
tion that arises from quark and gluon fluctuations. The mesonic resonances 
that dynamically form at lower energies, in particular in the s-channel of 
the four-fermi interaction, lead to effective low-energy degrees of freedom, 
the mesons. In the flow equation approach this is described by dynamical 
hadronization, see [291 EOl [221 EI]. This approach to QCD has been put 
forward for QCD in |12l I32j and is reviewed in [2]. 

The present dynamical PQM model is discussed in detail in [13], includ- 
ing its direct connection to the FRG approach to QCD. The flow equation 
for the free energy of QCD is shown in Fig. [H It is important to empha- 
size that the depicted flow is exact, no higher loop diagrams are missing. 
Here, Vk[4>\ stands for the effective action in the presence of an infrared 
cut-off scale k and the free energy is given by the effective action evalu- 
ated at the equation of motion. For each degree of freedom a fully dressed 
one-loop structure emerges. For example, the first two loops for the gluons 
and ghosts comprise the pure glue sector in QCD: the full gluon and ghost 
propagators are screened in the presence of dynamical quarks and mesons, 
respectively. 

An example for the matter contribution to the glue propagator is shown 
in Fig. [2] with the vacuum polarisation. The quark loop screens the gluon 
propagator, and in particular with increasing density the in- medium effects 
get larger. This has crucial effects, such as the reduction of the decon- 
finement transition temperature of the pure YM system when full QCD, 
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Fig. 2. Quark loop correction to the gluon propagator contributing to the matter 
back-couphng to the gauge sector present in full QCD. 



including the dynamical quarks and meson, is considered. 

The simple additive structure has important consequences for the PQM 
model. We can easily identify the different parts of the model in terms of full 
QCD fluctuations. For example, the fluctuating quark-meson sector of the 
model is given by the last two diagrams in Fig. [TJ In other words, removing 
the glue sector in Fig. [1] leaves us with a dynamical quark-meson model, 
for a review see [2jj. In turn, the glue and ghost loops in Fig. [H evaluated 
in the background of a Polyakov-loop give the Polyakov-loop potential 
in full QCD. This potential only agrees with the Yang-Mills potential if we 
neglect the matter fluctuations, such as the vacuum polarisation depicted 
in Fig. [21 This entails that with Fig. [2] we have access to the change of the 
Polyakov-loop potential in the presence of matter fluctuations. 

It is precisely this simple structure which has been used in [7] for a 
phenomenological estimate of the matter back-coupling to the glue sector. 
On the practical level, the back-coupling changes the dynamical Yang-Mills 
scale Aym, present in the gluon and ghost propagators, to the dynamical 
QCD scale Aqcd- In the Polyakov-loop potential this scale manifests it- 
self in the critical temperature Tq. Therefore, in full QCD this parameter 
receives a flavour and chemical potential dependence via the vacuum polar- 
isation, see Fig. [21 leading to To{Nf,fj.) |7|. By now, we also have access to 
the full QCD potential at vanishing chemical potential derived from solving 
the fully coupled QCD flows [2lll21[33j. This result confirms the phenomeno- 
logical estimate in [7], and will soon be extended to finite chemical potential, 
thus resolving the current ambiguity in the Polyakov-loop potential in the 
Polyakov-loop extended models. 

The discussion above allows us to put forward the flow equation of QCD 
where the gluonic degrees of freedom have been integrated out. This leads 
to a flow for the free energy which only involves the last two loops in Fig. [H 
and the first two loops lead to the Polyakov-loop glue potential Ogiue • They 
also lead to modifications of the matter interaction which is included in 
the initial conditions in the quark-meson sector that capture the correct 
vacuum physics. Ignoring the additional modification of the quark and 
meson dispersion this leads us finally to the flow equation [16[ [T3] for the 
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QCD free energy r^qcD.k = ^k, 

MA(.>.*,4) = ^ {^coth (I) + Acoth (f ) (1) 

AN Nf 

^ [1 - N,{T, CD, ^) - iV5(r, ^; ^)] 

where A; is the infrared cut-off scale. The quasi-particle energies are given 
by, i = q,TT, a, 



Ei = -^fc2 + mj , ml = h^"^ , ml = 20'^ , ml = 20'^ + icP'^Q'l , (2) 

with (p = (o", vf) and primes denote derivatives with respect to 0^. The 
Polyakov-loop enhanced quark/anti-quark occupation numbers read 

Nq{T, /i; = ^ ^ 3^e(^9-/^)/^ + 3«>e2(^9-M)/T + (A{e^-^)IT 

and A''g(r,//;$,^) = iVg(r, -/i; ^, ^>). Eq. ([I]) as it stands is stih a flow 
equation for the free energy of full QCD with the approximation of a classical 
dispersion of quarks and mesons. The gluons have been integrated out, 
their dynamics is directly stored in the glue Polyakov-loop potential rigiue 
as well as the initial conditions for Qa- Hence, only the matter part rimatter,fc 
depends on the cut-off scale k, to wit SfcriqcD.k = Sfcrimatter.fc- This leads 
us to the final expression for the QCD free energy Oqcd = ^QCD,fc=0) 

+^^matter,A(0-,7f,^>,^) . (4) 

The flow equation ([1]) is indeed that of the matter sector of the two-flavour 
QCD flow derived in [12] with <^> = ^{Aq). At sufficiently large initial cut-off 
scale A the initial matter part of the free energy, f^mattcr.A, is just a local 
Yukawa- type action of quarks and mesons, its parameters are fixed with 
phenomenology in the vacuum with /t^, and m^. We also remark that 
the independence of the full free energy from the initial cut-off scale A, that 
is Sa^^qcd = 0, enforces A dependent terms in r^matter.A- These terms can 
be determined from the flow at A, see the reviews [2 H [22 } [23 | [2 ^ [25] . 
Phenomenologically, they can be understood as the high-energy part of the 
vacuum fluctuations, see |34[ 135] . 

In the following we present results based on the full integration of the 
flow of r^QCD without any further approximation: the flow ([T]) is solved on 
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Fig. 3. Phase diagram with constant Tq (left) and with running To(/i) (right). 



a (T-grid for general fixed backgrounds Note in this context that the 

couplings in the initial condition r^QCD.A are insensitive to the choice of 

By mapping out the ^, ^-plane this provides us with the final result 
f^QCD(o-,7r, ^>,^) = OQCD,k=o(o",7f, ^>,^). This goes beyond the approxi- 
mation used in [13] where the flow was solved on the mean-field solutions 
<I>((t), ^((t). It has been argued there that the latter already provides a quan- 
titatively reliable approximation to the full solution and the present results 
fully confirm this argument. More details will be presented elsewhere |36| . 
The present computation allows us to solve the equations of motion for $, 
$ and a without further approximations and hence to determine the free 
energy and other thermodynamic quantities. 

The above discussion of the Polyakov loop potential only exemplifies how 
results obtained within the first principle FRG approach to QCD can be used 
in order to constrain model computations. It is by no means restricted to the 
current example and provides a way of systematically improving the models 
towards full QCD. Moreover, lattice results provide further input and in 
particular also help to improve the quantitative precision of the QCD flows. 
In our opinion, in combination this opens a systematic path of mapping out 
the phase diagram of QCD. 

3. Phase structure at physical pion masses 

Fig. [3] summarizes our findings for the two-flavour phase structure at 
physical pion masses, ttt-tt = 138 MeV, and for a constant Tq (left) and 
running Tq{^) parameter (right). The gray band denotes the width of the 
temperature derivative of the Polyakov-loop (solid black line) at 80% of its 
maximum height. The dashed black line labels the inflection point of the 
conjugate Polyakov-loop. Both Polyakov-loop transitions, in the following 
referred to as deconfinement transition, are crossovers at low /i. With a 
constant Tq (left panel), we find an almost ^ independent deconfinement 
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temperature which is shghtly below the chiral crossover (dotted blue line). 
At /i ~ 200 MeV the chiral and deconfinement transitions start to deviate 
and we find a region in the phase diagram where chiral symmetry is partially 
restored and confinement still persists. The picture changes when we include 
the matter back-reaction to the gluonic sector via a running To(/i) (right 
panel). In this case, both transitions lie close to each other throughout the 
whole phase diagram and no room for a chirally restored and confined phase 
is left. Remarkably, a similar scenario was found in another recent two- 
flavour non-perturbative functional study with Dyson-Schwinger equations 



With or without matter back-coupling to the YM system, we always 
find a critical endpoint around (T^^^, /i'^^P) ^ (20 . . . 30, 290) MeV. The 
endpoint is located at low temperatures, which results from the inclusion of 
quark and meson fiuctuations. In standard mean-field calculations, where 
the mesonic fiuctuations are ignored, T^^^ is typically much higher [38]. 
Already the inclusion of the fermion vacuum fiuctuations lead to much lower 
T'^^^ values, see e.g. [391 ES [M] • 

Note however, that the present approximation lacks accuracy at large 
chemical potential. There we expect baryonic degrees of freedom to be 
important. We indeed envisage that it is not so much the baryonic off- 
shell fiuctuations but rather their importance for the true ground state 
that matters. Diquark fiuctuations, however, may play a quantitative role. 
These considerations are supported by the findings in two-colour QCD, see 
e.g. [ini HU HJl 133]. Hence we conclude that the present approximation 
lacks predictive power for ^i/T ~ 2. In turn, for smaller ratios it is reliable 
and as our computations show, there is no sign of a critical point in this 
regime. More details on this matter will be presented in [36]. 




In the chiral limit, chiral symmetry is exactly restored at high temper- 
atures and/or densities and the chiral transition is a phase transition and 
not a crossover, see Fig. H] (left panel). At /Li = we encounter a second- 
order phase transition, in agreement with 0(4)-universality arguments |44j . 
At large chemical potential, we find an interesting phase separation: the 
second-order chiral transition splits into two branches for decreasing temper- 
atures. The inner branch (at smaller chemical potential) shows the inward- 
bending behaviour characteristic for studies including thermal and quantum 
fiuctuations. It turns into a first-order transition for smaller temperatures 
with a critical endpoint. The outer branch bends outwards and stays of 
second-order. This behaviour is in agreement with previous FRG findings 
in the pure quark- meson model in the chiral limit, cf. [39j . Interestingly, the 




4. Chiral limit 



8 



MaxBorn28 printed on February 6, 2012 




critical temperature of the endpoint T'^^^(/i^^^) is almost independent of 
the explicit symmetry breaking, i.e., independent of the pion mass. This is 
in contrast to T^^^(/i = 0) which for decreasing pion mass also decreases. 
Hence, the endpoint lies in the same temperature range as for physical pion 
masses. In turn, the critical chemical potential /i^^^ changes with decreas- 
ing pion mass to smaller values. This can be understood as follows: the 
chemical potential at which the phase transition line hits the /x-axis is re- 
lated to the value of the quark mass and also to the sigma meson mass. 
Both masses are smaller in the chiral limit, cf. also |45j . 

The transition splitting is accompanied by two minima in the sigma me- 
son mass as a function of the chemical potential. This is shown in Fig. [5] 
(right panel) for a temperature above the endpoint. In contrast to the 
sigma meson mass, the pions stay massless until chiral symmetry is com- 
pletely restored, which happens at the outer transition branch. For larger 
chemical potential, both meson masses degenerate due to the restored chiral 
symmetry. 



5. Small pion masses 

As previously shown, the chiral limit phase structure has peculiar fea- 
tures at large chemical potential that are not present at physical pion mass. 
In the following, we demonstrate how the chiral limit is connected to the 
physical mass point and study the influence of an increasing pion mass on 
the chiral transition splitting in the phase diagram. 

For a non-vanishing explicit chiral symmetry breaking parameter, the 
second-order chiral transition at moderate chemical potential turns into a 
crossover. Apart from this, the phase structure remains qualitatively the 
same. At small pion masses, for example m.,^ = 50 MeV, the splitting in the 
chiral transition persists and we find a critical endpoint on the inner branch. 
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Fig. 5. Similar to Fig. m but for = 50 MeV. 

again around Tc « 20 — 30 MeV, cf. Fig. [5] (left). For smaller temperatures 
the transition is of first-order. The outer branch remains a crossover for all 
temperatures. The sigma meson mass as a function of chemical potential 
shows a behaviour similar to the chiral limit. The second minimum is still 
visible, but much weaker, see Fig. [5] (right). For a larger explicit symmetry 
breaking this effect is washed out. Furthermore, the inflection point of the 
chiral condensate, which determines the transition line, is also smeared out 
at larger pion masses. At physical pion masses no splitting is observed 
anymore, cf. Fig. [3l 

6. Pressure 

In order to understand the properties of the chiral splitting region in the 
phase diagram in more detail, we also compute the pressure in this regime. 
In Fig. [B] we show the pressure normalized by its Stefan-Boltzmann value 
as a function of ^ for m-,^ = 50 MeV at the fixed temperature T = 30 
MeV. Clearly, the pressure in the splitting regime is not a monotonically 
rising function. A local minimum in the scaled pressure appears at about 
the chemical potential where the sigma meson mass has its local maximum. 
Only for chemical potential beyond the second transition branch the pres- 
sure increases significantly. Note however, that the dropping of the pres- 
sure for these densities is unphysical and hints at some shortcoming of the 
Polyakov-loop implementation in PQM/PNJL models. As already pointed 
out above, the coefficients of the effective Polyakov-loop potential are fitted 
to lattice data at vanishing chemical potential. Some density aspects of the 
matter back-coupling are included by our modifications where Tq has been 
replaced by T^^Nf, fj,), but a fully dynamical treatment within the FRG is 
still missing. 

Indeed, the slope of the used effective Polyakov-loop potential ansatz 
is known to be steeper than the corresponding potential in full QCD, see 
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Fig. 6. Pressure in the splitting region, T = 30 McV, for — 50 MeV. 

\12\ [53l [2] which might result in such an unphysical behaviour. In any case 
we observe that a small change in the parameters of the Polyakov loop 
potential lead to qualitative changes in the regime under discussion. Hence, 
in order to get a more realistic description of the QCD phase diagram, in 
particular in this density regime, the full dynamics of the glue potential 
Ogiuc has to be taken into account. In our opinion this, together with 
the inclusion of baryonic and diquark degrees of freedom already discussed 
above, will give us access to the highly interesting nuclear matter regime at 
low temperature and relatively large densities. 



In the present contribution we have discussed the effects of quantum and 
thermal fluctuations as well as the pion mass sensitivity of the phase struc- 
ture of QCD. This is done within the dynamical Polyakov-extended quark- 
meson model where fluctuations are included via the functional renormaliza- 
tion group. We have also argued that the PQM model can be understood as 
a well-controlled approximation to first principle QCD. The precise knowl- 
edge of the approximations involved can be used to systematically improve 
the PQM model towards full QCD by successively fixing the model param- 
eters with QCD input. Specifically, we have considered the back-coupling 
of quarks to the glue sector of QCD which results in a Nf- and /i dependent 
modification of the Tq parameter in the Polyakov-loop potential. 

At physical pion masses, the FRG phase diagram exhibits a critical 
endpoint whose location is pushed towards large chemical potential and low 
temperature in comparison with a standard mean-field approximation. This 
emphasizes the significant influence of fluctuations. We could rule out the 
emergence of an endpoint for small chemical potential below /i/T ~ 2. 

In the chiral limit, a splitting of the chiral transition line at large chem- 
ical potential and below T < 50 MeV could be confirmed within our trun- 



7. Conclusion and Outlook 
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cation. Two transition branches with different transition orders emerge for 
smah temperatures. On the inner branch, a second-order critical endpoint 
as the endpoint of the first-order transition fine could be determined. Within 
the splitting region chiral symmetry is still broken. 

Interestingly, the area of the splitting region is not very sensitive to the 
pion mass. For a small pion mass the splitting phenomenon still persists and 
the critical temperature of the endpoint is almost pion mass independent. 
For non- vanishing pion masses the chiral transition changes from the chiral 
0(4) universality to a crossover. Since fluctuations weaken the chiral phase 
transition the crossover is in addition washed out and disappears when we 
approach physical pion masses. However, some remnants of the second 
transition branch are still present in the chiral condensate at physical pion 
masses. 

The investigation of the thermodynamics in the splitting region might 
shed more light on the nature of the different transitions and the emergence 
of a critical endpoint. Concerning the pressure in this region we found 
that the Polyakov-loop potential has some deflciencies at low temperature 
and large chemical potential. Within the splitting region, the slope of the 
normalized pressure as a function of chemical potential becomes negative, 
which is an unphysical behaviour and hampers the computation of other 
thermodynamic quantities there. However, in this regime of the phase dia- 
gram baryonic degrees of freedom are of importance but have been neglected 
so far. 

Furthermore, the inclusion of the matter back-coupling to the glue sec- 
tor favours the coincidence of chiral and Polyakov-loop transition lines also 
at larger chemical potential. This is in clear contradistinction to the mean- 
field results where one observes a separation of the two transitions. This, 
together with large- A'c investigations, has been used as an indication for a 
quarkyonic phase, where chiral symmetry is restored but matter still con- 
fined. Given the intricacies related with the interpretation of the Polyakov 
loop as an order parameter for the confinement /deconfinement phase transi- 
tion at finite density, we merely note that our results with fiuctuations cast 
some doubt on the support provided by the mean-field investigations with 
Polyakov-extended models. Seemingly, it is the lack of matter back-coupling 
in these investigations which triggers the separation of the transition lines. 
Note also that this missing back-coupling is built-in in the large- iVc limit. 

In summary, the dynamical Polyakov-loop extended quark-meson model 
provides a good approximation to full QCD at small densities. For densities 
beyond ///T > 2 the present approximation has to be extended to include 
baryonic degrees of freedom, the effects of in-medium propagation and that 
of multi-scatterings in a dense medium. 
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